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ABSTRACT 

We present VLT-SINFONI K-band integral field spectroscopy of the central galaxies in the 
cool core clusters A1664, A2204 and PKS 0745-191, to probe the spatio-kinematic properties 
of the Paa and ro-vibrational H2 line emission. 

In A 1664 the two emission-line velocity systems seen in our previous Ha spectroscopy appear 
in both Paa and H2 emission, with notable morphological differences. The recession velocity 
of the red component of Paa increases linearly with decreasing radius, particularly along an 
8 kpc filament aligned with the major axis of the underlying galaxy and the cluster X-ray 
emission. These kinematics are modelled as gravitational free-fall as gas cools rapidly out 
of the hot phase. In A2204 the gas shows 3 or 4 filaments reaching radii of 10 kpc, three of 
which lie towards 'ghost bubbles' seen in X-ray imaging by Sanders et al. For PKS 0745- 
191, we confirm the twin-arm morphology in the narrow -band images of Donahue et al.; 
the Paa kinematics suggest rotational motion about an axis aligned with the kiloparsec-scale 
radio jet; on nucleus, we find an underlying broad Paa component (FWHM 1700 km ) and 
a secondary H2 velocity system redshifted by +500 km s^^ . 

The H2 v=l-0 S(3)/Paa ratio is highest in the most isolated and extended regions where 
it matches the levels in the NGC 1275 filaments as modelled by Ferland et al. Regions 
with much lower ratios highlight active star formation and are often kinematically quiescent 
(FWHM < 200 km s^^ ). Our findings suggest that the three clusters may be captured in dif- 
ferent stages of the 'cold feedback' cycle of Pizzolato & Soker, with A1664 in a short-lived 
phase of extreme cooling and star formation prior to an AGN heating event; PKS 0745-191 in 
an outburst state with the AGN accreting from a cool gas disk, and A2204 in a later phase in 
which cool gas is dragged out of the galaxy by the buoyant rise of old radio bubbles. 

Key words: galaxies:clusters:individual: A1664 - galaxies:clusters:individual: A2204 - 
galaxies:clusters:individual:PKS 0745-191 - cooling flows - intergalactic medium 



1 INTRODUCTION 

The past decade has witnessed a flood of new observational insights 
into the X-ray cooling flow phenomenon in the cores of galaxy 
clusters. On the X-ray side, results from XMM-Newton and Chan- 
dra have led to a sharp downward revision in X-ray cooling rates 
and exposed a strong deficit of line emission from gas cooling be- 
low one third of the ambient cluster temperature (see review by 
Peterson & Fabian 2006). The implication is that the cooling of the 
hot X-ray gas is quenched, with the most likely cause being that 
episodic outbursts from radio-loud active galactic nuclei (AGN) 
somehow heat the gas (see review of McNamara & Nulsen 2007). 



* Based on observations performed at the European Southern Observatory, 
Chile (Programme ID: 77.A-0057(A)). 



Significant progress has also been made on the discovery and 
characterization of various phases of cool gas and dust in cluster 
cores. In CO line emission. Edge (2001) reported detections of 16 
clusters, consistent with 10^ - lO"'^ Mq of H2 at 20-40 K for 
a standard C0:H2 conversion (see also Salome & Combes 2003). 
Interferometry showed further that the CO emission is localised 
within the central few arcsec of the central cluster galaxy (CCG) 
(Edge & Frayer 2003; Salome & Combes 2004). The CO emission 
often coincides with K-band rovibrational emisison from much 
smaller masses (~ 10^ - 10^ Mq) of H2 at ~ 2000 K (Edge 
et al. 2002; Jaffe, Bremer & van der Werf 2001). Similarly, mid- 
infrared spectroscopy with Spitzer has revealed pure rotational 
H2 emission from molecular gas at 300-400 K in several systems 
(e.g. Egami et al. 2006a, Johnstone et al. 2007). All these molec- 
ular emission line components appear to correlate well with the 
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strength of Ho? emission whose prevalence in these systems, often 
in the form of spectacular filamentary systems, has been known for 
several decades. The most comprehensive spectroscopic survey for 
Ha emission is that of Crawford et al. (1999), who showed that 
in the most luminous systems the emission is powered by massive 
star formation. The latter is also the heating mechanism for the ex- 
cess thermal far infrared dust emission discovered by Spitzer MIPS 
photometry of a sample of CCGs by Egami et al. (2006b). These 
results are corroborated by the more extensive Spitzer photometric 
survey of line-emitting CCGs by Quillen et al. (2008), revealing 
excess far-infrared emission in about half the sample, consistent 
with star-formation in the known molecular gas reservoirs (O'Dea 
et al. 2008). 

The emerging consensus is that whilst AGN heating success- 
fully offsets most of the cooling in a time-avergaed sense, residual 
cooling to sub-X-ray temperatures persists at levels of, at most, a 
few per cent of the rates previously expected. This gas accumu- 
lates in molecular gas reservoirs within the central galaxy, giving 
rise to star formation and associated emission line nebulosity. The 
physical processes at work in these cooled gas reservoirs are, how- 
ever, not well understood. It is not known whether the star forma- 
tion is steady and continuous or subject to outbursts, perhaps trig- 
gered by the central radio-loud AGN or interaction with a passing 
galaxy. The excitation mechanisms of the emission line nebulosity 
are also uncertain: whilst star formation can power the Ha emis- 
sion on energetic grounds, the overall emission line spectrum dif- 
fers markedly from a typical HII region or starburst galaxy, with 
optical forbidden and rovibrational H2 lines being far stronger than 
observed in such systems (e.g. Wilman et al. 2002). Simulations by 
Ferland et al. (2008,2009) suggest that the H2 emission in CCGs 
may arise from excitation by cosmic rays or heating via dissipative 
magnetohydrodynamic processes in dense gas shielded from exter- 
nal radiation fields. Support for this scenario comes from Hubble 
Space Telescope observations of the Ha filaments in NGC 1275 
by Fabian et al. (2008), who showed that the filaments can be sup- 
ported by threads of magnetic field and thereby stabilised against 
collapse and star formation. 

To shed further light on the cool gas systems in cluster cores, 
we have begun detailed studies of individual systems using optical 
and near-infrared integral field spectroscopy. The aim is to map the 
warm ionized gas and hot molecular hydrogen to reveal their spatial 
distributions, kinematics and excitation properties. In the first pa- 
per, Wilman, Edge & Johnstone (2005) presented high-resolution 
near-infrared integral field spectroscopy (IPS) of the nuclear re- 
gion of NGC 1275 using UIST on UKIRT, leading to the discov- 
ery of a 50-pc radius disk of hot H2 and a dynamical measurement 
of the nuclear black hole mass, the latter in good agreement with 
the Mbh — o" relation. The existence of this disk has been con- 
firmed by first results from the adaptive optics Near-infrared Inte- 
gral Field Spectrograph (NIFS) on the Gemini telescope (McGre- 
gor et al. 2007). Observations of NGC 1275 with the Smithsonian 
Millimetre Array by Lim et al. (2008) have shown how this disk is 
fed from the inflow of cool molecular gas along filaments extend- 
ing out to radii of 8 kpc. These filaments coincide with the coolest 
X-ray plasma between the two radio cavities, and together these ob- 
servations demonstrate, for the first time, how the cooling of X-ray 
gas in a cluster core feeds the central supermassive black hole. 

In a second study, Wilman, Edge & Swinbank (2006) pre- 
sented VLT-VIMOS optical IFS of four of the most Ha-luminous 
CCGs known. The full 2-d view revealed a variety of disturbed 
morphologies, ranging from the smooth but distorted to the ex- 
tremely clumpy and disturbed, with velocity gradients and split- 



tings of several hundred kms~^ on scales of 20 kpc or more, and 
evidence for an association between disturbed Ha emission and 
secondary galaxies. The global Ha kinematics match those of the 
CO(l-O) emission in single-dish data, suggesting that they arise in 
the same gas clouds. The kinematic coupling between the Ha and 
ro-vibrational H2 emission appears less strong. Despite the dis- 
turbed kinematics, the ratio of optical forbidden line emission to 
Ha is remarkably constant within and between CCGs, implying 
that the line ratios in these extreme environments are saturating due 
to widespread star formation. Less luminous systems exhibit more 
spatial variation in line ratios (e.g. Hatch et al. 2007), as expected 
if star formation is more patchy and irregular. 

To build on this work, we present in this paper near-infrared 
IFS with VLT-SINFONI of three CCGs (two from the Wilman et 
al. 2006 sample), to perform a direct comparison between the Paa 
and H2 emission at high spatial and spectral resolution. Compared 
with our previous studies, VLT-SINFONI instrument offers a 5-fold 
improvement in spatial sampling over VLT-VIMOS, and substan- 
tially better sensitivity than UKIRT- UIST 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of the CCGs of A 1664, A2204 and PKS 0745- 
191 were taken in service mode in 2006 with the SINFONI inte- 
gral field unit (IFU) (Eisenhauer et al. 2003) on UT4 of the Very 
Large Telescope (VLT) at the European Southern Observatory in 
Paranal, Chile. The observation dates and exposure times are listed 
in Table 1. The IFU was operated throughout in the 8x8"^ 
field-of-view mode, with 32 slitlets of 0.25" width, with a detec- 
tor pixel scale of 0.125" along the slits. The K-band grism span- 
ning a wavelength range of 1.95-2.45/im at resolution ~ 4000 was 
employed, resulting in a dispersion of 2.45 x W^'^pm per pixel. 
The instrument was operated without adaptive optics and with an 
ABBA pattern of 300 s on/off source exposures to enable sky sub- 
traction. The seeing was approximately 1" throughout. Standard 
pipeline data reduction was carried out using the ESOREX pipeline. 
In brief, the reduction comprises the following steps: correction 
for bad pixels, flat fielding, correction for geometrical distortion, 
arc-line-based wavelength calibration, subtraction of sky emission 
(using ABBA sequences) and merging of the individual frames 
into a single data cube for each object, with sub-pixel dithering 
enabling re-sampling to square pixels of 0.125". A refined back- 
ground subtraction was performed during post-processing in IDL, 
using a background spectrum synthesised from regions free of CCG 
line emission in the pipeline-reduced datacubes. All subsequent sci- 
entific analysis was also performed in IDL. The spectra were flux 
calibrated but a telluric standard was not observed, potentially cor- 
rupting the spectra outside the 2.08-2.34/^m window of clean at- 
mospheric transmission. This is for the most part not a problem 
because most of the emission lines of interest lie in this range. The 
exception is PKS 0745-191, where Paa and H2 v=l-0 S(5) fall just 
blue-ward of 2.08/im. In this one target, we corrected for this dur- 
ing the analysis by dividing the raw spectra by a scaled version of 
the standard Mauna Kea atmospheric transmission curve for a com- 
parable spectral resolution (the data were obtained from the UKIRT 
worldwide webpages, and originally produced using the program 
IRTRANS4). The scaling (i.e. the absorption optical depth) was 
adjusted in order to flatten out the continuum across the absorption 
troughs. 

In the analysis we devote one section to each of the three 
CCGs, followed by some generalised conclusions. The cosmo- 
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Table 1. Observation log 



Target 


Redshift 


Observation 


On-source exposure 


Scale 






dates 


time (s) 


(kpc arcsec"^) 


A 1664 


0.1276 


2006 Apr 4, May 19,20,30,31 


7200 


2.3 


A2204 


0.1514 


2006 Jul 1,4,5,7, Aug 2,3 


10800 


2.6 


PKS 0745-191 


0.1028 


2006 Apr 5,6,20,26,27 


10200 


1.9 



Mpc"\ Q.M = 0.3 and 



logical parameters Hq — 70 km s 
SIa =0.7 are assumed throughout. 



3 RESULTS ON A1664 

A1664 has a redshift of z = 0.1276 and from the Spitzer mid- 
infrared spectral energy distribution O'Dea et al. (2008) inferred a 
star formation rate of 14.6 Mq yr~^ , compared with a cooling rate 
deduced from X-ray spectra of 60^20 yr~^ ■ At ~ 50 kpc, the 
size of the star formation region implied by the Schmidt law is the 
second largest in the O'Dea et al. sample (most of the rest are in 
the 10-20 kpc range). However, the source is poorly resolved in the 
mid-infrared which would instead imply a size around 15 kpc. 

A 1664 was the most spectacular and complex of the sample 
of four CCGs in our earlier Ha study with VLT-VIMOS (Wilman, 
Edge & Swinbank 2006). The Ha emission was characterised by a 
31 kpc filamentary structure passing through the nucleus in a NE- 
SW direction, an orientation shared by the larger-scale cluster X- 
ray emission. Two velocity components were detected along the 
full length of the filament with a splitting of almost 600 kms~^ in 
the nuclear regions. A bipolar nuclear morphology in the blue kine- 
matic component suggested that some kind of starburst-driven out- 
flow may be responsible for this complexity. However, the presence 
of a small galaxy just beyond the end of the filament (and within 
lOOkms"^ in line-of-sight velocity) led us to favour a scenario 
in which this companion cluster galaxy had plunged through the 
cluster core, disturbed the gas reservoir and triggered a burst of star 
formation. The current VLT-SINFONI IFU data enable us to map 
this structure at higher spatial resolution through the Paa emission, 
which is an isolated emission line and thus free of the line blending 
seen in the Ha-l-[NII] complex. 

The Paa line (rest-wavelength 1.8756/im) in each 0.125" cell 
was fitted with a gaussian emission line complex atop a flat con- 
tinuum over the wavelength range 2.07505-2. 13630pim (i.e. 250 
spectral bins). Chi-squared values were computed for fits with zero, 

O 

one or two velocity components: Xo^ Xi X2- For Xo ~ Xi > 30 
one component was fitted; for Xi "Xi > 10, two components were 
deemed necessary. A similar procedure was used to fit the H2 v=l-0 
S(3) emission line (rest- wavelength 1.9576^m). 



3.1 Gas morphologies and kinematics 

An immediate insight into the complex kinematics is shown in the 
velocity distributions derived from the emission line fits to Paa 
and H2 v=l-0 S(3) (Fig. Q}. Two distinct velocity systems are ev- 
ident: one spanning -400 to -1-200 km s^^ , the other ^00 to - 
600kms~^ (relative to a systemic redshift of z — 0.1276). Here- 
after we refer to these as the 'red' and 'blue' components respec- 
tively. The morphologies of these components are shown alongside 
a K-band continuum map in Fig.|2l 



Velocities of cells with fitted Pa-aipiia in A1664 (wrt z=0.1276) 
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Figure 1. Velocity distributions for Paa and H2 v=l-0 S(3) in A1664 
as a function of spatial pixel number within the IFU aperture, 
clearly demonstrating the existence of two distinct velocity sys- 
tems: a red system (denoted by the crosses) and a blue system 
(diamonds). The dividing line between the two is chosen to be - 
400 km s"^ . 



As expected, the Paa total intensity map bears a close resem- 
blance to the Ha map shown in Wilman et al. (2006). Minor dif- 
ferences exist between the morphologies of the 'blue' and 'red' 
velocity components in Paa and Ha, but these are due to the im- 
proved separation which is now possible for Paa compared with 
the Ha-l-[NII] complex. The reduced extinction in the near-infrared 
compared to the optical may also lead to morphological differences 
(Crawford et al. (1999) derive reddening of £(5 — 1^) — 0.46 mag 
from the B aimer decrement). 

The morphologies in Paa and H2 v=l-0 S(3) are notably dif- 
ferent, as was already hinted at by the disparity in their line widths 
in our original K-band long-slit spectrum (Edge et al. 2002). The 
peaks of all the emission line maps are spatially offset from one 
another and from the K-band continuum nucleus of the galaxy. In 
Fig. ID we show an HST WFPC2 (F606W filter) image of A 1664 
(private communication, CP. O'Dea), overlaid with contours of the 
Paa and H2 v=l-0 S(3) total intensities. This wide V-band image 
reveals a complex morphology with several bright knots of star for- 
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mation and obscuring dust lanes. The Paa emission peaks on the 
central dust lane and for the most part follows the distribution of 
the star-forming knots. The H2 emission largely avoids the star- 
forming region to the south of this dust lane, but exhibits a sec- 
ondary peak coincident with another prominent dust feature. A full 
analysis of the HST image will appear elsewhere. Spectra for a 
selection of spatial regions are shown in Fig. |4] and confirm the 
marked variations in H2 v=l-0 S(3)/Paa ratio, to which we will 
return in section 6. 

The two-dimensional gas kinematics are shown in Fig.|5] The 
most coherent large-scale motions are those of the Paa red com- 
ponent, whose recession velocity increases as one moves in to- 
wards the nucleus. The strongest velocity gradient is along the 
north-east-south- west direction ('axis 1'), where the FWHM of 
the emission is also elevated with respect to the other parts, at 
~ 400 kms^^ . The blue Paa component is largely quiescent, with 
FWHM < 200 km s"^ in the regions with the most coherent ve- 
locity structure. The H2 kinematic maps are necessarily noisier but 
the red velocity field resembles that of Paa. 

Integrating the emission within the SINFONI field-of-view, 
the total line luminosities are 5.3 x 10*^ erg s"^ (Paa) and 3.1 x 
10*^ ergs^^ (H2 v=l-0 S(3)). These are 5 and 12 times larger, re- 
spectively, than the line luminosities implied by the fluxes quoted 
in Edge et al. (2002) within a 1.22" -wide north-south slit. 



3.2 Interpretation: inflowing filaments of cooling gas? 

The increasing recession velocity of the red component of Paa with 
decreasing projected radius could be interpreted as the acceler- 
ated infall of gas from the outer regions of the galaxy towards 
the nuclear regions. For more insight into this. Fig. |6] shows one- 
dimensional cuts in this velocity field along three separate axes in- 
dicated in Fig.|5] These demonstrate that, beyond an initial starting 
radius, the velocity increases linearly with decreasing radius. For 
axis 1, the velocity increases by 500 km s^^ over 5 kpc, along axis 
2 by 250 km ^ over 3kpc, and along axis 3 by 400 km ^ over 
4 kpc. 

Similar kinematics were observed Lim et al. (2008) for the 
CO filaments in NGC 1275 extending out to 8 kpc, and in the near- 
infrared long-slit spectroscopy of Wilman et al. (2000) for Cygnus 
A. Lim et al. demonstrated that the linear increase of velocity with 
radius could be modelled as free fall within the gravitational poten- 
tial of the galaxy, for a parcel of cool gas released at rest from some 
initial radius. The CO filaments in NGC 1275 are co-spatial with 
the coolest X-ray gas in the system (orientated roughly perpendic- 
ular to the radio lobes), suggesting that the filaments represent gas 
which has cooled out of the hot phase, decoupled from hydrostatic 
equilibrium, and started to free fall inwards. 

The soft X-ray emission in A 1664 (Fig. 9 of Wilman et 
al. 2006; Kirkpatrick et al. 2009) is also aligned with the inflow 
of gas along the Ha filament (axis 1), and with the major axis 
of the K-band continuum. The alignment of the major axis of the 
CCG with the cluster X-ray emission in A 1664 was noted by Allen 
et al. (1995). Following Lim et al. (2008), we model the velocity 
field by considering free-fall in a gravitational potential given by 
the Hernquist (1990) profile characteristic of elliptical galaxies: 
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where Ad is the total mass of the galaxy and a a scale radius. 
As shown by Lim et al., the radial velocity v{r) of a test particle 
released from rest at radius r = ro is given by: 



The parameters Al and a are not well-constrained for A 1664, 
so we adopt fiducial values of M = lO^'^ M© and a — 13.3 kpc. 
The latter assumes an effective radius Re = 24 kpc (the mean of a 
broad distribution measured by Schneider et al. (1983) for a CCG 
sample) and the relation Rc — 1.8153a given by Lim et al. Anal- 
ysis of the K-band IFU image (Fig.|2ll and the HST image (Fig. [5] 
ignoring the disturbed central regions), loosely constraints the ma- 
jor axis iic to the range 18-27 kpc, consistent with the assumed 
value. 

For starting radii ro = 12, 8 and 4 kpc. Fig. |7] shows the re- 
sulting 1-d profiles of line-of-sight velocity versus projected radius 
when the filament is oriented at various angles, 9, to the sky plane. 
The steepening of the model velocity gradient as the filament is 
oriented closer to the line-of-sight implies that the observed ve- 
locity profiles for axes 1,2 and 3 cannot be modelled by varying 
the viewing angle for a fixed starting radius. Within the context of 
this model, we find that these gradients can be accounted for with 
the following combinations of parameters (ro, 6*): (12kpc, 65°) for 
axis 1; (4kpc, 40°) for axis 2; (8kpc, 55°) for axis 3. These are 
not to be regarded as best fit parameter values, as there are clearly 
degeneracies between the parameters (ro, 0) in addition to the un- 
certainties due to the poorly constrained values of R^ and M. 

What is not readily explained by this model, however, is the 
high observed FWHM for the infalling gas along axis 1 . This may 
suggest that the filament is observed from a more foreshortened 
perspective with the infalling gas covering a range in solid angle. It 
may instead arise from turbulence due to hydrodynamic interaction 
of the clouds with the diffuse medium. 

In the absence of better constraints on the parameters of the 
gravitational potential, this analysis serves only to demonstrate that 
an infall interpretation of the observed velocity field is plausible. 
The observed kinematics may have a different origin, e.g. from 
an interaction with a secondary cluster galaxy, as suggested by 
Wilman et al. (2006). A second alternative is that the filamentary 
Ha emission and extended 'bar-like' X-ray emission (Kirkpatrick 
et al. 2009) may be due to a 'sloshing' motion of the CCG within 
the cluster potential, as invoked to explain the X-ray and Ha fil- 
ament in A1795 (Fabian et al. 2001). In that case, however, the 
Ha filament has quiescent kinematics over the bulk of its length 
with more violent motions only seen around the powerful radio 
source within the CCG (Crawford et al. 2005). 



4 RESULTS ON A2204 

In Wilman et al. (2006), A2204 was found to have an irregular 
Ha morphology consisting of a several weakly-defined filaments 
extending from the nuclear regions, possibly coincident with dust 
features seen in an Hubble Space Telescope optical image. The 
velocity field appeared to comprise 3 clumps of emission at ve- 
locities of -200, -50 and -100 kms"^ (relative to z = 0.1514), 
with a FWHM~ 200 km s~ , rising to 800 km s~ ^ in the vicinity 
of the nucleus. From the VIMOS data, the secondary galaxy lo- 
cated ~ 10 kpc in projection to the south-west of the CCG was 
measured to have a line-of-sight velocity of ~ — lOOkms"^ rel- 
ative to the CCG. It is thus plausible that the two galaxies are in- 
teracting. Narrow-band Ha-l-[NII] imaging was presented by Jaffe 
et al. (2005), which revealed diffuse line emission out to radii of 
24 kpc with a north-west-south-east alignment. 

O'Dea et al. (2008) derived a mid-infrared estimate of the star 
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Figure 2. Continuum and emission line maps for A1664. The maps for H2 1-0 S(3) have been smoothed with a gaussian of FWHM 0.5". For each of the 
two emission lines, the red, blue and total intensity maps are plotted on a common greyscale range for ease of comparison. The peak surface brightness of 
Paa maps is 4 times that of the smoothed H2 v=l-0 S(3) maps. The crosshair in each panel denotes the location of the continuum nucleus of the galaxy, as 
assessed from the K-band image. 



formation rate of 14.7 M© yr~ . The latest X-ray and radio obser- 
vations are consistent with a cooling rate of 65 M© yr^'^ down to 
the lowest detectable temperature (Sanders, Fabian & Taylor 2008). 
Within 20 kpc of the X-ray cluster centroid, Sanders et al. iden- 
tify between 5 and 7 X-ray surface brightness depressions, each 
~ 4 kpc in size. None of them is associated with the radio emis- 
sion, so they may be 'frothy ghost bubbles'. Their estimates of 



the enthalpy contained within the bubbles fall short by an order 
of magnitude of the amount required to offset cooling in the central 
100 kpc. 

The present IFS data delineate the gas filaments much more 
clearly, showing that they extend roughly perpendicularly from the 
nuclear regions to the north, south, east, with an additional protru- 
sion to the west-north-west. The latter feature is absent in Ho v=l- 
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Figure 3. An HST WFPC2 image of A1664 (F606W filter) overlaid with contours of the emission in Paa (left) and H2 v=l-0 S(3) (right). 
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Figure 4. Spectra of selected regions in A1664 showing tlie spatial variations in Paa and H2 line ratios. The y-axis units on all spectra are 

10~^^ erg cm^^ s^^fim^^ . 



S(3) emission, but the overall morphologies in Paa and H2 v=l- 
S(3) are otherwise very similar. Intriguingly, three of these ex- 
tended features appear to lie along the radial vectors towards 'ghost 
bubbles' identified by Sanders et al. The bouyant rise of these bub- 
bles may have dragged the cool line-emitting gas out of the central 



regions, as seen most spectacularly in the Perseus cluster (Fabian 
et al. 2003). 

The kinematic maps in Fig. |9] reveal that the extended fila- 
ments are blue-shifted in velocity from the emission in the main 
body of the galaxy. This suggests that they may have been pulled 
out of the galaxy towards the observer; the velocity gradient ap- 
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Figure 5. Kinematic maps for the red and blue velocity components of Paa and H2 v=l-0 S(3) in A1664. All velocities are in kms ^ , with reference to a 
systemic redshift of z = 0.1276. For reference, the crosshair locates the continuum nucleus. 



pears to change sign at the tip of the northern filament, suggest- 
ing that the outward motion stalls there. Such velocity reversals 
have been observed in the Ho filaments in NGC 1275 by Hatch et 
al. (2006) and are in qualitative agreement with the expected flow 
patterns around rising buoyant bubbles. The line-width of Paa rises 
to 1000 km FWHM on nucleus but there is no evidence for 
the presence of distinct broad and narrow-line components (cf. 
PKS 0745-191). 

Spectra are shown for various regions in Fig.[TO] As in A1664, 
significant variations in H^fPaa line ratio are seen: in extended 
'filamentary' regions A, B and E, H2 v=l-0 S(3) is stronger than 
Paa, whereas in region F H2 v=l-0 S(3) is extremely weak. The 
spatially-integrated line luminosities are 3.0 x 10*^ ergs^^ (Paa) 
and 2.8 x lO'*^ ergs"^ (H2 v=l-0 S(3)); these are 2.2 and 3.7 times 
larger, respectively, than those implied by the slit fluxes quoted in 
Edge et al. (2002). 



5 RESULTS ON PKS 0745-191 

Donahue et al. (2000) presented Hubble Space Telescope imaging 
of the central cluster galaxy of PKS 0745-191 in emission lines of 
H2 v=l-0 S(3) and Ha-l-[NII]. Our Paa and H2 v=l-0 S(3) emission 
line maps shown in Fig. [TT] are in excellent correspondence with 
their observations. The extended line emission is characterised by 



two roughly-perpendicular, ~ 5 kpc-long, arm-like features which 
protrude on the western side of the nucleus, and by a short spur 
which extends ~ 1.5 kpc from the nucleus on the eastern side. As 
in A2204, the Paa emission is more strongly nucleated than the 
H2 emission. The current estimates of the star formation and X- 
ray cooling rates are 17.2 M0 yr~^ and 200^3^ Mq yr"^ , respec- 
tively (O'Dea et al. 2008). 

The H-band continuum imaging of Donahue et al. revealed 
a putative jet-like feature extending about 1" from the nucleus in 
a westerly direction, which they observed to be coincident with a 
similar feature in the 2cm radio image of Baum & O'Dea (1991). 
The feature is also marginally present in our K-band image in 
Fig. [TT] The 6cm radio image of Baum & O'Dea revealed ex- 
tended emission but no clear indication of any classical jet or 
lobe features. The negative residuals (after substraction of regu- 
lar isophotes) in the H-band imaging revealed a dust morphology 
consisting of several 'fingers' 1-2" in length emanating from the 
nucleus; the strongest of these coincides with the spur of line emis- 
sion to the east-south-east of the nucleus. 

Kinematic maps are shown in Fig. [12] They reveal that the 
two arms of emission exhibit distinct velocities approximately 
100 km ^ apart. There is evidence for a rotational motion, partic- 
ularly in the Paa velocity field, about an east-west axis well-aligned 
with the position angle of the emerging radio jet. The rotation curve 
derived from a reconstructed 0.25" wide pseudo-slit shows a sharp 
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Figure 6. Velocity cuts for the red velocity component of Paa in A 1664 as a function of projected nuclear distance along the three axes shown 
in Fig. [5] Negative nuclear distance corresponds to the numbered ends of the axes shown in Fig. [5] To give some indication of the errors on 
these velocities, the different symbols refer to three adjacent cells at the same nuclear distance along each axis. 




Figure 7. For 3 values of the starting radius, the solid lines show the expected radial free-fall velocities for a test particle released from rest in a model 
gravitational potential, as described in section 3.2. The dashed lines in each panel show, from right to left, the corresponding Une-of-sight velocities which 
would be observed as a function oi projected radius when these filaments are oriented at angles of 25, 35, 45, 55, 65, 75, and 85 degrees to the plane of the sky. 



jump in velocity across the nucleus, despite the modest (1") see- 
ing. This may well be a rotating disk seen close to edge-on, analo- 
gous to that observed in NGC 1275 by Wilman et al. (2005); future 
studies under better seeing and with adaptive optics could use it to 
constrain the mass of the black hole. The connection between the 
two arms of extended emission and this putative disk is not clear: 
the two arms could represent cooling material falling into this disk, 
although there is no clear velocity gradient along the arms (any mo- 
tion could be in the plane of the sky, or drag forces may have caused 
inf ailing clouds to reach a terminal velocity); alternatively, the arms 
may represent ambient material pushed aside by the emerging radio 
source. The spectra for various regions shown in Fig.|13|reveal that 
the H2/Paa line ratio exhibits some spatial variation, but not to the 
same extent as in A 1664 and A2204. 



The spatially-integrated line luminosities are 
4.2 X 10"^ erg s"^ (Paa) and 2.0 x 10"*^ erg s"^ (H2 v=l- 
S(3)); these are 5.3 and 6 times larger, respectively, than those 
implied by the slit fluxes quoted in Edge et al. (2002). 



5.1 Evidence for a second velocity component: a dusty 
jet-driven outflow? 

The velocity dispersion maps for Paa and H2 v=l-0 S(3) both show 
peaks within 0.5" radius of the nucleus. Examination of the binned 
spectra in Fig.ll3lshows that on nucleus (region D) the H2 lines pos- 
sess a red-shoulder which compromises the fit with a single gaus- 
sian. Closer inspection of the spectra of the four central 0. 125" pix- 
els (see Fig.ll4t shows that this is due to a distinct H2 emission line 
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Figure 8. Continuum and emission line maps for A2204. For ease of comparison, the crosshair in each panel denotes the location of the continuum nucleus 
of the galaxy, as assessed from the K-band image, and the emission line maps use a common greyscale. The galaxy to the south-west is a secondary cluster 
galaxy with a line-of-sight velocity blue-shifted by 100 km ^ relative to the central cluster galaxy (see Wilman et al. 2006 for details). The dashed lines on 
the 1-0 S(3) image show the radial vectors in the directions of the 7 'ghost bubbles' identified in X-ray imaging by Sanders et al. (2008) (6 of the bubbles lie 
at projected radii of 20kpc, outside the SINFONI field-of-view). 
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Figure 9. Kinematic maps for Paa and H2 v=l-0 S(3) in A2204. All velocities are in kms ^ , with reference to a systemic redshift of ^ = 0.1514. For 
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system offset in velocity by ~ 500 km s^^ . More tentatively, a 
second component is also apparent in the spectrum of region E at 
a relative redshift of ~ 330 kms~^ . It is most pronounced in the 
para-H2 lines of v=l-0 S(2) and v=l-0 S(4), whilst the ortho-lines 
of v=l-0 S(l), S(3) and S(5) have more symmetric profiles. Phe- 
nomenologically, this could arise if the apparent ortho:para ratio in 
the secondary component were suppressed significantly below the 
value of 3 which pertains for local thermodynamic equilibrium at 
T > 200 K. This could occur under UV fluorescent excitation (e.g. 
Sternberg & Neufeld 1999) but in this case we would expect the 
higher order v=2-l H2 lines produced in the radiative cascade to be 
stronger relative to the v=l-0 lines. However, this test is inconclu- 
sive since the peak in v=2-l S(3) coincides with the lower velocity 
component (with v=2-I S(3)/v=l-0 S(l) ~ 0.2) and any redshifted 
component lies below the noise. More plausibly, and consistent 
with the interpretation of the strong H2 filament emission in CCGs 
(see section 6), the second velocity component could be produced 



in dense, well-shielded gas at 10-100 K in which the H2 emission is 
excited by collisions with secondary non-thermal particles, as mod- 
elled by Ferland et al. (2008,2009) under the 'cosmic ray' case. 
In such cool gas, the LTE H2 ortho:para ratio would be substan- 
tially below 3. The v=2-I/v=l-0 line ratios would be comparable to 
the UV fluorescence case, as set by the balance between radiative 
cascade and coUisional de-excitation rates (the latter have recently 
been revised substantially, as discussed by Ferland et al. 2009). 
Regardless of the excitation mechanisms, these secondary veloc- 
ity components may be related to some form of outflow triggered 
by the emergence of the kiloparsec scale jet visible in the 2cm radio 
map. 

A similar asymmetry is not apparent in the Paa profiles, but a 
closer inspection of the profiles of this line in the nuclear regions 
reveals the presence of an underlying broad component centred at 
the same redshift as the narrow component (Fig. |15ll. The broad 
component has a FWHM of 1700 km . Such a line width falls 



Integral field spectroscopy of ionized and molecular gas in cool cluster cores 1 1 



PKS 0745-191 K-band continuum 



PKS 0745-191 Pa-alpha intensity 



PKS 0745-191 H2 1-0 S(3) intensity 



i 



10 kpc 



12 3 4 



5 

Arcsec 



6 7 8 9 1 2 3 



4 5 
Arcsec 



7 8 9 




Figure 11. Continuum and emission line maps for PKS0745-191. For ease of comparison, the crosshair in each panel denotes the location of the continuum 
nucleus of the galaxy, as assessed from the K-band image, and the emission line maps use a common greyscale. 
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Figure 12. Kinematic maps for Pao and H2 v=l-0 S(3) in PKS 0745-191. All velocities are in km s^^ , with reference to a systemic redshift of z = 0.1028. 
For reference, the crosshair locates the continuum nucleus. Also shown is a 1-d velocity cut for a 2-pixel (0.25") wide pseudo-slit aligned north-south through 
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Figure 13. Spectra of selected regions in PKS 0745-191. The y-axis units on all spectra are 10 erg cm ^^m ^. 



in the middle of the distribution derived by Landt et al. (2008) from 
their near-infrared spectral survey of a broad-line AGN sample. The 
inclusion of the additional broad component clearly improves the 
overall fit to the emission line profile, but does not significantly 
affect the FWHM and velocity fields derived from the single com- 
ponent fit in Fig. 1121 



6 EXCITATION OF THE GAS 

Table 2 lists the H2 v=l-0 S(3)/Paa emission line flux ratio for 
the labelled spatial regions in each galaxy. This particular line ratio 
was chosen because both lines are within the observed wavelength 
range for all three targets and are sufficiently close together to min- 
imise uncertainties arising from differential flux-calibration errors 
and extinction. 



As already hinted by cruder long-slit spectroscopy (e.g. Edge 
et al. 2002), the nuclear regions correspond to local minima in 
the H2/PaQ ratio, most likely due to enhanced Paa emission from 
a photoionized AGN component. The regions with the highest 
H2/Paa ratio are generally isolated features morphologically dis- 
tinct from the contigious body of emission within the galaxy, viz. 
the extended filaments in A2204 (regions A, B and E); the extended 
'arms' in PKS 0745-191; the extremity of the infalling gas stream 
in A1664 (region E). The H2 v=l-0 S(3)/PaQ ratios in these regions 
are extremely high, typically within 30 per cent of unity. The latest 
modelling by Ferland et al. (2008,2009) accounts for these high val- 
ues (and other line ratios which are anomalous with respect to pho- 
toionization predictions) by invoking extra heating from magneto- 
hydrodynamical waves and/or cosmic-ray-induced non-thermal ex- 
citation of the molecular lines. The models invoke emission from 
gas at a range of densities pinned to magnetic fields at a common 
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Figure 14. The top four panels show spectra of the four nuclear pixels in PKS 0745-191 (within the cross-hair in region D of Fig.[T3]l. The 
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S(4) at an offset of ~ 330 kms~^ (dotted lines). The expected positions of the two velocity components in the H2 v=l-0 S(l), S(3), S(5) and 
v=2-l S(3) lines are indicated by the dot-dashed lines. 



pressure (P/k — nT) matching that of the diffuse X-ray plasma. 
The emergent spectrum depends on the assumed power-law density 
distribution and can be arbitrarily tuned to match any one line ra- 
tio, but for the 'Horseshoe filament' in NGC 1275 the best-fitting 
Ferland et al. (2009) models for the 'extra heat' and 'cosmic ray' 
yields H2 v=l-0 S(3)/Paa values of 1.37 and 0.73, respectively. 
These values bracket the measured ratios in the most extended re- 
gions of our galaxies, suggesting that the same physical processes 
are at work. 

There are, however, several spatial regions with extremely low 
H2 v=l-0 S(3)/Paa ratios in the range 0.15-0.3, namely regions 
A, B, C (red velocity component), D (blue) in A1664, region F 
in A2204, and region F in PKS 0745-191. Star formation is most 
likely the dominant excitation mechanism in these regions. Indeed, 
region A in A 1664 coincides clearly with a knot of star formation 
emission in the HST image of Fig.jS] To further illustrate the role of 
star formation, we have made a comparison with the near-infrared 
spectroscopic survey of a sample of ultraluminous infrared galax- 
ies (ULIRGs) by Dannerbauer et al. (2005). This unbiased sample 
of 24 ULIRGS was drawn from those observed with the Infrared 
Space Observatory (ISO), and the vast majority have mid-infrared 



and optical spectral classifications of starburst or HII region. Us- 
ing the published emission line fluxes in Table 2 of Dannerbauer 
et al., we compute average H2 v=l-0 S(3)/Paa ratios of 0.19 (me- 
dian), 0.3 (mean) with a standard deviation 0.46 for the 22 spec- 
tral regions with detectable H2 v=l-0 S(3) and Paa (some of the 
ULIRGs are split into independent spatial regions). A further 8 of 
their spectra exhibit Paa but no detectable H2 v=l-0 S(3) and in 
some other cases the latter line may blend with [SiVI] 1.957; the 
average H2 v=l-0 S(3)/Paa ratios quoted above are thus likely to 
over-estimate the true value, but are nevertheless comparable with 
the lowest values measured in our CCGs. 

Regions with H2 v=l-0 S(3)/Paa intermediate between the 
above extremes are likely to reflect a blend of star formation and 
non-photoionization 'filament-excitation'. A crude insight into the 
physical location of the star formation within the CCG can be ob- 
tained by substracting 1//fil times the H2 v=l-0 S(3) image from 
the Paa image, such that the residual Paa emission mainly re- 
flects the contribution from star formation. The value of /fil is 
set by the maximum observed value of H2 v=l-0 S(3)/Paa for 
each galaxy, which we take to be characteristic of the 'filament- 
excitation', i.e. /fil = 1.3 for A2204, 0.75 for A1664 and 0.86 
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Figure 15. Profiles of the Paa emission line for the central four pix- 
els along the north-south pseudo-slit in PKS 0745-191 (from which 
the rotation curve of Fig.[T2]is derived). Fits with single and double 
gaussians are overlaid, showing the existence of a broad-line com- 
ponent. The individual panels are labelled with their distances from 
the nucleus in 0.125" pixels along the pseudo-slit. 



for PKS 0745-191. The residual Paa intensity maps are shown in 
Fig. [16] For A2204, the star formation is confined to the main body 
of the galaxy with an additional clump 4kpc to the north-west; 
for A 1664, there appear to be five main clumps of star-formation: 
two circunmuclear components ~ 1 kpc either side of nucleus (one 
in the red velocity component, one in the blue), a second compo- 
nent ~ 3.7 kpc south of the nucleus in the blue velocity compo- 
nent, region A (see Fig. |4j, and an area ~ 3.5 kpc to the north 
north-west of the nucleus. These features mostly correspond with 
the star-forming regions visible in the HST image (Fig. [3]l, except 
for the two circumnuclear regions which are obscured by the dust 
lane. The infalling gas stream (along axis 1 in Fig. |5j is relatively 
devoid of star formation. With reference to the kinematic maps, 
most of these regions correspond to gas with the lowest velocity 
dispersion {FWHM < 200 kms"^). For PKS 0745-191, the 
residual image emphasises the main body of the galaxy and the 
area between the 'arms' of H2 emission. These residuals are not 
a pure star formation tracer, since some of the nuclear emission 
with depressed H2 v=l-0 S(3)/Paa may arise from AGN-related 
photoionization. When the previously-quoted total Paa luminosi- 
ties of the three CCGs are converted to Ha luminosities assuming 
case B recombination (Ha/Paa= 8.45), and then to a star forma- 
tion rate (SFR) using the Ha-SFR relation of Kennicutt (1998) 
[SFR(MQyr~^) = L(Ha)/1.26 x iC^ergs-M, we obtain the 
following SFRs: 36MQyr~^ (A1664), 2OM0yr"^ (A2204) 
and 28MQyr~^ (PKS 0745-191). These values exceed the in- 
frared measurements of O'Dea et al. (2008) of 14.6, 14.7 and 
17.2MQyr^^, respectively. If we correct for the contribution 
of filament excitation to the Paa emission (assuming H2 v=l- 
S(3)/Paa= 0.15 for star formation and the above values of /fil 
for each CCG), the implied SFRs become: 11 (A1664), 6 (A2204) 
and 15 Mq yr"^ (PKS 0745-191). A correction for dust extinction 
of j4v = 2 — 6 mag towards the Paa-emitting star-forming re- 
gions (the range estimated from a local sample of luminous infrared 
galaxies by Alonso-Herrero et al. 2006) would boost the SFRs by 
30-100 per cent, and the infrared and Paa estimates would then be 
more comparable. 

These substantial variations in H2 v=l-0 S(3)/Paa contrast 
with the constancy of the ratio of optical forbidden line emission 



to Ha in A 1664 and A2204, as reported by Wilman et al. (2006). 
In the above scenario where the observed line emission reflects a 
spatially-varying blend of photoionization due to star formation 
(and possibly an AGN) and non-photoionization 'filament' exci- 
tation, this implies that the optical line ratios for the two types of 
emission are very similar, whilst the infrared H2 v=l-0 S(3)/Paa ra- 
tios are manifestly quite different. In both cases, the optical forbid- 
den lines and the bulk of the HI lines are produced in gas in a rela- 
tively narrow range of parameter space with T ~ 10* K and density 
n ~ lO'^'^ cm^'^ . In the 'cosmic ray' case of filament excitation 
favoured by Ferland et al. (2009), the volume emissivities of the 
HI and H2 lines extend deep into the much cooler, denser and neu- 
tral parts of the clouds, with a secondary peak in Ha emissivity at 
100-1000 K due to coUisional excitation by fast electrons and a sec- 
ondary peak in ro- vibrational H2 emission at 10-100 K. However, 
the volume occupied by this dense gas in their model is much lower 
and it contributes negligibly to the observed optical line emission, 
which is dominated by the partially ionized gas. Moreover, Ferland 
et al. show that the internal extinction across a 0.3 pc thread (sev- 
eral hundred of which compose each 70 pc-wide filament in NGC 
1275) is A\r ~ 17 mag, arising mostly within the dense molecular 
regions. This implies that the observed optical emission is domi- 
nated by the low-extinction partially ionized gas. In contrast, the 
K-band extinction is Ak ~ 0.1 Ay, enabling the Paa and H2 lines 
to probe much deeper; the broad emissivity distributions of these 
lines renders the observed H2 v=l-0 S(3)/Paa ratio sensitive to the 
details of the gas density distribution, leading to spatial variations 
in H2 v=l-0 S(3)/Paa over and above those caused by a varying 
mix of 'filament' and star formation (photoionization) emission. 



7 DISCUSSION 

Recent analyses by Rafferty et al. (2008) and Cavagnolo et 
al. (2008) have demonstrated the existence of a threshold for the 
onset of star formation and radio-loud AGN activity in CCGs when 
the central X-ray cooling time drops below 5 x 10* yr (equivalent 
to a central entropy threshold of 30keV cm^ ). This underscores 
the role of AGN heating as the ultimate regulator of X-ray cooling 
in these environments. The existence of a common (sharp) thresh- 
old in cooling time for star formation and AGN activity suggests 
that both are a direct consequence of X-ray cooling. In the light 
of this, we discuss the current VLT-SINFONI observations and the 
prospects for IFS studies of larger samples of CCGs. 

In A 1 664, the evidence for the gravitational free-fall of cooled 
Paa and H2-emitting gas from radii ~ 10 kpc suggests that the sys- 
tem is observed at an early stage in this particular cooling/heating 
feedback cycle. Freshly cooled gas rains down on the CCG and 
gives rise to extensive star formation within the molecular cloud 
lifetime of ~ 10^ yr (Mouschovias et al. 2006), which is compa- 
rable to the free-fall timescale from 10 kpc. The AGN has either 
not been fuelled yet or its activity is still insufficient to arrest the 
X-ray cooling. Some of the infalling gas will eventually fuel the 
AGN directly if its angular momentum is low, otherwise accretion 
may proceed via a disk (e.g. as observed in NGC 1275; Wilman et 
al. 2005). The central radio source is indeed weak with the unre- 
solved NVSS 1.4 GHz detection of 36mJy implying a monochro- 
matic power of lyLiy ~ 2 x 10'"' ergs^^ , placing it at the lower 
end of the radio power distribution for CCGs below the cooling 
time (entropy) threshold in Cavagnolo et al. These observations of 
free-falling cool gas and a weak AGN are consistent with the 'cold 
feedback' model of Pizzolato & Soker (2005). In this model, blobs 
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Figure 16. Residual Paa images formed by subtracting a scaled version of the H2 v=l-0 S(3) image from that of Paa. The residual image 
should emphasise the regions where the Paa emission is due to star formation, as described in section 6. 



Table 2. Spatially-resolved emission line ratios within the CCGs (f indicates the nuclear region). 



Galaxy 


Region 


H2 v=l-0 S(3)/Paa 


A1664 


A 


0.28 ± 0.025 




B (blue) 


0.16 ±0.017 




B (red) 


0.19 ± 0.08 




C (blue) 


0.75 ± 0.042 




C (red) 


0.17 ± 0.086 




D t (blue) 


0.23 ± 0.026 




D t (red) 


0.41 ± 0.020 




E 


0.72 ± 0.057 


A2204 


A 


1.29 ±0.09 




B 


1.23 ±0.075 




ct 


0.69 ± 0.018 




D 


0.81 ± 0.056 




E 


1.25 ± 0.14 




F 


0.32 ± 0.054 


PKS 0745-191 


A 


0.86 ± 0.037 




B 


0.57 ± 0.039 




C 


0.62 ± 0.028 




Dt 


0.36 ± 0.01 




E 


0.47 ± 0.015 




F 


0.31 ±0.017 
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of gas overdense by at least a critical factor 5c — 2 relative to 
the ambient medium initially cool slowly and flow steadily inwards 
over a period of ^ 10* yr since the last AGN heating event. Over 
the next ~ 10^ yr they enter a free-fall plunge into the CCG over 
the final ~ lOkpc of radius as they cool catastrophically, contract 
in size and drag forces become negligible. Some of this cool gas 
may form stars and some may fuel the AGN, leading to the next 
major heating event and the repetition of the cycle. A1664 may 
well be captured in this final phase with an another AGN heating 
event due within ~ 10^ yr. A Chandra study of A1664 by Kirk- 
patrick et al. (2009) confirms that the system is in a low-state of 
AGN activity insufficient to offset the X-ray cooling, which oc- 
curs at a rate comparable to the star formation rate. The alignment 
between the elongated X-ray bar (possibly a collapsed cavity sys- 
tem; Kirkpatrick et al.) and the Ha filament further strengthens the 
hypothesis that the star formation is a direct result of the X-ray 
cooling and not triggered by interaction with a secondary galaxy as 
suggested by Wilman et al. (2006). 

Not surprisingly, neither of the other two CCGs in our sam- 
ple is in this short-lived pre-outburst state with cool gas free-falling 
onto the CCG. With its kiloparsec-scale radio jet and higher power 
(vL^ ~ 9 X lO^ergs"^ from the NVSS 1.4 GHz detection), 
PKS 0745-191 is arguably observed shortly after an AGN outburst 
with the central black hole still being fed from the gas disk ob- 
served in the extended Paa emission; on nucleus, the broad-line 
Paa component and secondary H2 velocity system point to active 
gas accretion and a possible outflow. A2204, with its 'ghost cavi- 
ties' and ~ 10 kpc-scale radio source (Sanders et al. 2008), is prob- 
ably observed at a somewhat later post-outburst epoch. As argued 
by Rafferty et al., star formation within the cool gas reservoir can 
persist for more than 10* yr and exhibit detectable signatures over 
several feedback cycles and will thus not appear synchronized with 
the latter. Future IFU observations of a larger sample could identify 
more of these 'pre-outburst' systems, by for example targeting the 
subset of clusters in the Cavagnolo et al. sample which have central 
cooling times below the critical threshold but weak (i.e. currently 
undetected) radio emission. Studies of the multiphase cool gas in 
such systems could shed light on the black hole fuelling mecha- 
nisms and further constrain the 'cold feedback' model. 
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